Astronomy & Astrophysics manuscript no. nieppola'astroph 


February 5, 2008 


(DOI: will be inserted by hand later) 





Spectral Energy Distributions of a Large Sample of BL Lacertae 

Objects 

E. Nieppola* 1 , M. Tornikoski 2 and E. Valtaoja 1 ' 3 



1 Tuorla Observatory, Vaisalantie 20, FIN-21500 Piikkio, Finland 

2 Metsahovi Radio Observatory, Metsahovintie 114, FIN-02540 Kylmala, Finland 

3 Dept. of Physical Sciences, University of Turku, FIN-20100 Turku, Finland 



Received / Accepted 

Abstract. We have collected a large amount of multifrequency data for the objects in the Metsahovi Radio Observatory BL 
Lacertae sample and computed their spectral energy distributions (SED) in the log v - log vF - representation. This is the first 
time the SEDs of BL Lacs have been studied with a sample of over 300 objects. The synchrotron components of the SEDs 
were fitted with a parabolic function to determine the synchrotron peak frequency, v peak . We checked the dependence between 
luminosities at several frequency bands and synchrotron peak frequency to test the blazar sequence scenario, which states that 
the source luminosity depends on the location of the synchrotron peak. We also calculated broad band spectral indices and 
plotted them against each other and v pea k. 

The range of v peak in our study was considerably extended compared to previous studies. There were 22 objects for which 
\ogv pea k >19. The data shows that at 5 GHz, 37 GHz and 5500 A there is negative correlation between luminosity and v pea k, 
whereas in X-rays the correlation turns slightly positive. There is no significant correlation between source luminosity at syn- 
chrotron peak and v pea t- Several low radio luminosity-low energy peaked BL Lacs were found. The negative correlation between 
broad band spectral indices and v peak is also significant, although there is substantial scatter. Therefore we find that neither a TX 
nor a ro can be used to determine the synchrotron peak of BL Lacs. On the grounds of our results we conclude that the blazar 
sequence scenario is not valid. In all our results the BL Lac population is continuous with no hint of the bimodality of the first 
BL Lac samples. 

Key words, galaxies: active - BL Lacertae objects: general - Radiation mechanisms: non-thermal 



1. Introduction 

BL Lacertae (BL Lacs) are a subclass of active galactic nu- 
clei (AGN). They are characterized by the lack of strong emis- 
sion lines, rapid variability at all wavelenghts and strong polar- 
ization. Their spectral energy distribution (SED), in the log v 
- logvF -representation, consists of a synchrotron compo- 
nent at lower frequencies and an inverse Compton component 
at higher frequencies. The peculiar traits of the BL Lac class 
are most likely caused by Doppler-boosted radiation emanat- 
ing from a relativistic jet aligned close to the line of sight 
JUrrv & Padovanill995h . 

Traditionally BL Lacs have been discovered in either ra- 
dio or X-ray band, which led to their classification as radio- 
selected (RBL) and X-ray-selected (XB L) BL Lacs. B est- 
known RBL samples include the Uy iStickel etalJll99lt) . S4 
lStickel&KuhiHl994 and S5 faihr & Schmidtlll990h sam- 
ples an d among the mo st important XBL samples are the 
EMSS JGioia et alJll990l:lstocke et alJll99ll) and Slew Survey 
JPerlman et alll99fj|) . 

* e-mail: eni@kurp.hut.fi 



The two classes have different properties: RBLs are 
more variable, more luminous at radio and optical wave- 
lengths and have a higher polarization JStocke et alJ Il985t 
Llannuzi e t al J 19941 ). XBL s have a higher starlight fraction, 30- 
50 % (IStocke et al.lll985l) . and their morph ology is less core- 
domin ated in the radio than that of RBLs JPerlman & Stockel 
1993). Due to these differences they were initially regarded as 
separate classes of AGN. However, in recent years samples in- 
cluding intermediate objects have been found in surveys which 
combine X-ray and radio obs ervations. These include the RGB 
(RASS-Green Bank) sample dLaurent-Muehleisen et alJl99^) 
the Deep X-Ray Radio B lazar Survey (DXRB S) JPerlman et al 
| l998[lLandt et al.l200lh and the REX Survey dCaccianiga eTal 
1999). Their discovery has strengthened the view that the BL 
Lac population is continuous and RBLs and XBLs represent 
the opposite ends of the continuum. The reason for such a con- 
tinuity lies in the cutoff frequency of the synchrotron compo- 
nent in the SED. The synchrotron peak frequency of RBLs 
is in the radio/IR ban d and for XBLs the peak is mostly 
in the UV/X-ray band dGiommi et aljfl995l) . The intermedi- 
ate BL Lacs (IBL) have their synchrotron peak in the opti- 
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cal wavelengths. This explains why they were not observed 
in the first surveys. Following the synchrotron cutoff model 
the terminology was brought to more physical ground: the 
RBL/XBL -division was replaced by division into low-energy- 
peaked BL Lacs (LBL) and hig h-energy-peaked BL Lacs 
(HBL) JPadovani & Giommill995E . Most RBLs are LBLs and 
most XBLs are HBLs, but not all. The class boundaries can be 
loosely defined as v peak ~ 10 13 ~ 14 Hz for LBLs, v pea k ~ 10 15 ~ 16 
Hz for IBLs and w « 10 1718 Hz for HBLs. 

iFossati et ail lll998h linked the shape of the SED and 
the synchrotron peak frequency to the source luminosity: the 
lower the peak frequency, the more luminous the source. This 
would mean that LBLs are intrinsically more luminous than 
HBLs. This sequencing is based on the absence of high- 
luminosity HB Ls and low-luminosity LBLs. However, recently 
iGiommi et alJ d2005t) reported the possible discovery of high- 
luminosity HBLs in the Sedentary Surve y and evidence of low- 
power LBLs has also been discovered JPadovani et alJl2003l 
ICaccianiga & Marchall2004l) . These findings are at odds with 
the trend presented b y Fossati et al. 

lOhisellinil ll 19991) suggested that there is a class of BL 
Lacs whose synchrotron peak lies at even higher frequencies 
than that of conventional HBLs, v pea k > 10 19 Hz. These ob- 
jects can be called ultra-high-energy synchrotron peak BL Lacs 
(UHBLs) (Giommi et alJl200ll) . Following the dependency of 
the SED shape and luminosity, UHBLs are thought to be ex- 
tremely faint at radio wavelengths which is why they have es- 
caped notice. Extensive y-ray observations are needed to un- 
ambiguously identify them. 

In this paper we take a new approach to studying the prop- 
erties of the BL Lac population. Our goal is to plot the SEDs of 
the Metsahovi Radio Observatory BL Lac sample, which com- 
prises nearly 400 objects, including objects from all the best- 
known surveys at radio and X-ray wavelengths. This way we 
can examine the population properties of a sample with a wide 
range of attributes instead of focusing on one or two limited 
surveys. 

The aim of this paper was to test both the continuity of the 
BL Lac population and the blazar sequence scenario, and to as- 
sign a SED-based classification to those objects that previously 
had none. A large database of flux measurements is collected 
and SEDs are plotted for more than 300 BL Lacs. Each object is 
classified as LBL, IBL or HBL. The relationship of synchrotron 
peak frequency and luminosities at several frequencies is also 
tested along with the properties of broad band spectral indices. 
Throughout this paper we assume Hq = 65kms~ 1 Mpc~' and 
Qo = 1. All statistical tests were carried out using Unistat 5.5 
software. 



2. The sample 

The Metsahovi BL Lac sample includes 381 objects se- 
lected from the Veron-Cet ty & Veron BL Lac Catalogue 
JVeron-Cettv & VeroiJEoOoT) . hereafter VCV2000, and 17 ob- 
jects from literature. Given the northern location of the 
Metsahovi observatory, the source with the lowest declination 
in the sample is PKS 2223-114 at 5 =-11:13:41. The list of 



sample sources (table [3j is published electronically 1 . A large 
part of the objects in VCV2000 are from well-known and well- 
defined samples such as the Uy, S4, S5, EMSS and Einstein 
Slew Survey. Also BL Lacs from the first release of DXRBS 
identifications (Perlman et al J 199a) are included. No selection 
criteria (other than declination) in addition to the ones in the 
original surveys were imposed on the sample. The aim was to 
examine the behaviour of an extensive sample containing all 
known BL Lacs up to the year 2000. 

According to VCV2000 the Metsahovi BL Lac sample can 
be further classified as follows: 63 % are confirmed, 3 % prob- 
able and 8 % possible BL Lacs, 14 % are objects with high op- 
tical polarization and 12 % (including the BL Lacs taken from 
literature) lack any subclassification. In the sample there are 6 
sources that are not i n the later editions of Veron-Cetty & Veron 
BL Lac Catalogues JVeron-Cettv & VeroiJEoOlll2003l) . These 
objects have been excluded from the data analysis performed 
in this paper. 



3. The data 

In order to plot as accurate SEDs as possible, a large amount of 
data from several wavelengths was collected. Because simul- 
taneous multifrequency flux measurements are not available, 
datapoints from different epochs were searched from databases 
and literature. 

3.1. Radio data 

The starting point in collecting radio data were the Metsahovi 
observations at 37 GHz from late 2001 to January 2004. The 
full data set and a more detailed analysis about the 37 GHz 
behaviour of the sample sources will be published in a forth- 
coming paper (Nieppola et al., in preparation for the A&A). Of 
the BL Lac sample 137 objects were detected at <x > 4. There 
were 255 BL Lacs that were not detected and 6 that had not 
been observed yet during the time mentioned. The limiting flux 
of the Metsahovi radio telescope is about 0.2 Jy under optimal 
weather conditions. 

Our group has also obtained flux density data at higher ra- 
dio frequencies from our observations with the Swedish-ESO 
Sub millimetre Telescope (SEST) between 1987 and 2003 at 
3 mm and 1.3 mm. For some objects we got data from the 
RATAN observatory at frequencies 2.3, 4.8, 7.7, 11.2, 21.8 
and 30 GHz (Tornikoski et al., in preparation for the A&A). 
Additional low frequ ency datapoints we re found in VCV2000 
and WGA-catalogue llWhiteetalJl996l) . 

The large majority of radio data were obtained from the 
Astrophysical Catalogues Support System (CATS) maintained 
by Special Astrophysical Observatory, Russia 2 . The search re- 
sults from CATS contained data from over 140 different cata- 
logues and more than 100 radio frequencies. 



1 www.edpsciences.org 

2 http://cats.sao.ru/. 
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3.2. Data from other frequency bands 

The IR datapoints are from CATS, originating from IRAS- and 
2MASS-catalogues. The wavelengths used are 1.25, 1.65, 2.0, 
2.17, 12, 25, 60 and 100 ^m. 

The optical data also are mainly from CATS. Some data- 
points from V-band were added from Don ato et alJ ll200ll) . 

The X-ray data are from Einstein- and ROSAT-catalogues. 
The majority is from WGACAT JWhite et alJ 1 19961) and 
RBSC-catalogue dVoees et alJ [l999). D atap oints w e re also 
collected from the following papers: Donato et al. | d200lh . 
Lamer et all dl996l) and lLaurent-M uehleise n et alJ (11999). 
EMSS-d ata were included as well JGioiaetalJ [l990; 
IStockeet alll99ll) . 

All of the data fr om gamma region ca me from the Third 
EGRET Catalogue jHartman et alJ Il999h . Such data were 
available for only 14 objects. 



4. Computing the spectral energy distributions 

For all sources with a sufficient number of datapoints a spectral 
energy distribution was plotted in the form log v-log vF. All 
frequencies used are observed frequencies, they have not been 
reduced to rest frame frequencies. The synchrotron component 
of the SED was fitted with a parabolic function 



y 



= Ax 1 + Bx + C 



(1) 



in order to determine the synchrotron peak frequency Vp ea k = 
-B/2A. The fitting was successful for 304 objects; the rest were 
too sparsely sampled. The decision whether or not to include X- 
ray datapoints in the fit was based solely on a visual estimate 
for each individual object. All the SEDs (fig.ll3> are published 
electronically. The synchrotron peak frequencies are shown in 
table|3]col. 5. 

We note that using a simple parabolic function in the fit- 
ting produces some error especially among HBLs. In their case 
the X-ray datapoints are typically included in the rising syn- 
chrotron component, and therefore the parabola peaks after the 
X-ray domain. In reality, the synchrotron peak is expected to 
occur at or very close to the X-ray wavelengths and the de- 
cline to be more rapid. Thus the peak frequencies of the most 
extreme objects can be exaggerated. 

The objects were assigned a LBL/IBL/HBL classifica- 
tion according to y„ m t. O n the basis of the criteria of 
IPadovani & Giommil 1 1995), we decided to draw the bound- 
aries as follows: for LBLs, log v pea k < 14.5, for IBLs 14.5 < 
log Vpeak < 16.5 and for HBLs log v pea k > 16.5. Thus the peak 
frequencies of LBLs stretch up to the optical region, IBLs peak 
in the optical and UV-bands and HBLs from soft X-rays up- 
wards. 

4.1. The distribution ofv peak 

When the objects were classified as LBL/IBL/HBL according 
to their v pea i as described in the previous section, the divi- 
sion resulted in the three classes being almost equal in size. 
There were 98 LBLs, 96 IBLs and 110 HBLs. The distribu- 
tion is smooth and decreases steadily towards the higher peak 




peak 



Fig. 1. The distribution of the synchrotron peak frequencies in 
the Metsahovi sample. 

Table 1. Division of observational BL Lac classes to physical 
ones. 





LBL 


IBL 


HBL 


RBL 


84% 


10% 


6% 


IBL 


22 % 


36% 


42% 


XBL 


21 % 


25 % 


54% 



frequencies (fig. [Q. This can be a real effect, suggesting that 
sources in which electrons are accelerated to extreme energies 
are intrinsically rare, or due to selection effects in surveys. The 
frequency interval v pea k = 10 13 ~ 14 Hz is the most populated. 

Table ^presents how the classification with respect to the 
observational band relates to the one based on the SED. Here 
RBL classification has been assigned to objects in the Uy 
and S4 surveys, IB L to objects in RGB and 200 mJy BL Lac 
jBondi et all200ll) surveys and XBL to objects in the Einstein 
Slew Survey. T his RBL/XBL classification has been adopted 
from lGiommi et all i 19951) . Among the 304 sources for which 
the SED could be plotted there were 31 RBLs, 115 IBLs and 
48 XBLs. Several objects got multiple classifications. 

Table^clearly demonstrates how surveys in the X-ray en- 
ergies are more prone to find low-energy BL Lacs than radio 
surveys are to find HBLs. A large majority of RBLs really are 
LBLs. In fact, the 6 % of RBLs that turned out to be HBLs are 
the radio-luminous Mrk 42 1 and Mrk 501. Objects in the RGB 
and 200 mJy samples are more likely to be HBLs than LBLs, 
and only a third of them are truly intermediate. Of XBLs only 
half are HBLs and over 20 % are actually low-energy BL Lacs. 
These figures are certainly affected by the arbitrariness of the 
dividing boundaries between the classes, but the overall trend 
is expected to remain. It results from the fact that the X-ray 
luminosities of the samples in question are roughly the same 
whereas radio luminosities differ greatly. 

4.2. UHBL candidates 

As fig.[2shows, there were several objects in the sample whose 
synchrotron peak frequency was extremely high. Usually ob- 
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Fig. 2. Radio luminosity at 37 GHz plotted against synchrotron Fig. 3. Radio luminosity at 5 GHz plotted against synchrotron 
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peak frequency. 



jects with \ogv pea k > 17 are considered as extreme, here the 
number of such objects was 80, approximately 26 % of all the 
fits. For 22 objects \ogv pea k > 19 (table [2} and for 9 objects 
even log v pec ^ > 20, corresponding to a peak energy of ~ 0.4 
MeV. The SEDs of these sources are mainly very sparsely sam- 
pled, typically with datapoints from radio, optical and X-ray 
bands, and should be treated with caution. With that said, we 
note that even for objects with logVp ea * >20 the datapoints 
fit very well on the rising parabolic function. We note again 
that the actual position of the peak is probably exaggerated by 
the use of a parabolic fitting function as mentioned in sect. 4. 
Therefore the peak frequencies of these objects cannot be con- 
sidered as definite. Peaking near the MeV-region these sources 
would be excellent candidates for y-observations. 

5. Correlation between v pea k and luminosity 

5.1. High radio frequency (37 GHz) 

Out of the three BL Lac classes, LBLs had 37 GHz detections 
for 81 % of the sources, IBLs for 36 % and HBLs for 12 %. 
This seems to indicate that most HBLs have a radio flux well 
below the flux limit of the Metsahovi telescope. This prompted 
us to study the possible correlation between the synchrotron 
peak frequency and 37 GHz source luminosity more closely. 

To calculate the luminosities, redshift data wa s 
collected from VC V2 000, IStocke et alJ dl99ll) 
lLaurent-Muehleisen et ail ll 1999b . iDonato et all ll200ll) 
and SIMBAD database 3 . For sources which had no redshift 
available we used z = 0.4. Both a detection at 37 GHz and 
synchrotron peak frequency v pea j were available for 132 
objects. 

When the 37 GHz luminosity is plotted against the syn- 
chrotron peak frequency (fig. [5}, the correlation is easily seen. 
LBLs are at the high-luminosity end of the plot and the lu- 
minosity drops towards the high-energy regime. According to 
Spearman Rank Correlation Test, there is significant negative 
correlation at the 99 % confidence level. 



http ://simbad. u- strasbg.fr. 



No real evidence of a population of low-luminosity LBLs 
or high-luminosity HBLs is found. However, there are two 
LBLs with almost as low luminosities as HBLs disrupting the 
declining trend. In addition, the lowest luminosities are not 
claimed by objects in the high end of the peak frequency range, 
but by two IBLs with logv pea k < 16. The radio luminosities 
are widely scattered. An object with log vL = 44 can have peak 
frequency values ranging from \ogv pea j< = 13 to logv pea k = 17, 
approximately. Thus the radio luminosity cannot be used to de- 
termine the peak frequency of the source. 

5.2. Low radio frequency (5 GHz) 

The small number of HBL datapoints at 37 GHz convinced us 
to test the correlation also at 5 GHz. The number of available 
datapoints rose to 280. The correlation plot changed drastically 
(fig- B). The most noticeable difference is the appearance of 
several low-luminosity LBLs (lower left of the figure). They 
even reach lower luminosities than any of the HBLs. The over- 
all negative correlation is still present and significan t at 99 % 
level. This differs from the result obtained by IPadovani et alJ 
However, their DXRBS BL Lac sample consisted of 
only 31 objects mainly representing the LBL/IBL end of the 
plot. This is the region with most scatter in our figure and 
within this limited area the correlation is less obvious even with 
a larger number of datapoints. Only when the whole range of 
values of log v pea j is considered the trend becomes evident. 

When compar ed with the corresponding figure of 
iFossati et alJ Jl998l) our figure has much more scatter. A part 
of it is caused by a larger number of datapoints, but on closer 
examination our log Vnen k values for the Slew Survey are much 
more widely spread. In IFossati et alJ Jl998l) the Slew Survey 
objects have log v pea k — 15-19, while in our study they have 
\ogv pea k « 13-21. The Uy-sample takes on similar values of 
peak frequency in both studies. 

We note that there is no evidence of very high-luminosity 
HBLs. In fact, the extreme HBLs adopt quite intermediate 
luminosity values avoiding also the low-luminosity region. 
However, the spectra of possible high-luminosity HBLs would 
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Table 2. Objects for which log v peak >19. 



Source 


RA. (J2000) 


Dec. (J2000) 


log v peak 


1ES 0229+200 


02:32:48.6 


+20:17:17 


19.45 


RXS J03 14.3+0620 


03:14:23.9 


+06:19:57 


19.57 


2E 0323+0214 


03:26:13.9 


+02:25:14 


19.87 


2E 0414+0057 


04:16:52.4 


+01:05:24 


20.71 


1ES 0502+675 


05:07:56.1 


+67:37:24 


19.18 


EXO 0706.1+5913 


07:10:30.1 


+59:08:21 


21.05 


RXS J0847.2+1133 


08:47:12.9 


+ 11:33:52 


19.13 


1ES 0927+500 


09:30:37.5 


+49:50:25 


21.13 


RXS J 1008. 1+4705 


10:08:11.3 


+47:05:20 


19.67 


RXS J1012.7+4229 


10:12:44.3 


+42:29:57 


20.97 


EXO 1449.9+2455 


11:49:30.3 


+24:39:27 


19.83 


r(j 1218+304 


12:21:21.9 


+30:10:37 


19.14 


RXS J1319. 5+1405 


13:19:31.7 


+ 14:05:34 


20.85 


RXS J1341. 0+3959 


13:41:05 


+39:59:45 


20.06 


RXS J1353.4+5601 


13:53:28 


+56:00:55 


19.23 


RXS J1410.5+6100 


14:10:31.7 


+61:00:10 


20.25 


2E 1415+2557 


14:17:56.6 


+25:43:25 


19.24 


RXS J1456.0+5048 


14:56:03.7 


+50:48:25 


19.94 


RXS J1458.4+4832 


14:58:28 


+48:32:40 


21.46 


1ES 1533+535 


15:35:00.8 


+53:20:37 


19.68 


RXS J1756.2+5522 


17:56:15.9 


+55:22:18 


19.90 


RXS J2304.6+3705 


23:04:36.6 


+37:05:08 


21.01 
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Fig. 4. Optical luminosity at 5500 A plotted against syn- Fig. 5. X-ray luminosity at 1 keV and ROSAT band plotted 
chrotron peak frequency. against synchrotron peak frequency. 



very likely be totally featureless because of the powerful nu- 
cleus and thus the object would lack a redshift estimation. 
Therefore the redshift value assigned to featureless sources 
could have a big impact on the luminosity correlations and the 
appearance of high-luminosity HBLs (see section 5.6). 

5.3. Optical region (5500 A) 

Figure 0] shows the optical luminosity at wavelength 5500 A 
plotted against log v pea £. There is a significant, slightly negative 
correlation present at 95 % level. We see again that LBLs have 
more scatter than HBLs. Overall, the correlation is much less 
evident than in the case of radio luminosity. 



5.4. X-ray region 

In fig.[5]we have plotted the correlation between X-ray lumi- 
nosity and log Vpeak. In calculating the luminosities we have 
used both 1 keV data and data from the ROSAT band (0.1- 
2.4 keV). The error produced by the bandwidth difference is 
negligible when only the statistical properties of a large sample 
are considered. We note that in the case of X-ray luminosity, 
the correlation is positive and significant at 99 % level. This 
is in c ontrast with the blazar sequence scenario. iFossati et alJ 
(1998) state that the overall luminosity of HBLs is lower than 
that of LBLs (see their fig. 12). While they admit that in the 
X-ray band objects exhibit complex behaviour, the systematic 
rising trend presented by our findings is not predicted. 
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Fig. 6. Luminosity at the synchrotron peak frequency plotted Fig. 7. The distribution of log (S X /S r )- 
against synchrotron peak frequency. 



5.5. Peak luminosity L peak 

In addition to luminosities at defined frequency bands, we cal- 
culated for each source its luminosity at the synchrotron peak 
frequency. This is plotted against log v pea k in fig. [6] There is 
no significant correlation. Therefore we can decisively say that 
the source luminosi ty does not de pend on the synchrotron peak 
frequency. Fig. 7 of lFossati et alJ 09981) also shows the depen- 
dence of Lpeak and log v pea k- In their study there is a significant, 
yet weak correlation. Again, we remark that a larger number 
of datapoints and a wider range of log v pea k reveal the true be- 
haviour of the population and the lack of correlation. If only 
objects with log v pea k < 17 are considered, there is a weak neg- 
ative correlation also in our sample. On the other hand, when 
the high-energy tail with log v pea k > 17 is tested, we find a sig- 
nificant positive correlation. Therefore the distribution almost 
seems to assume a concave shape. 

5.6. The effect of z on luminosity correlations 

iGiommi et alJ d2005l) found numerous candidates for high- 
luminosity HBLs in the Sedentary Survey. All these objects 
seemed to reside at high redshifts (z >0.7). In our luminosity 
calculations we assumed z = 0.4 for featureless objects, which 
is a low value compared to Giommi et al. Using a too low red- 
shift value for a significant part of the sources would lead to 
a serious underestimation of luminosities. To take this into ac- 
count, we tested the luminosity correlations also by assigning 
redshifts z = 0.7, z = 1 and z = 1 .5 to those objects that had none. 

The effect was most pronounced in the case of radio lu- 
minosities at 5 and 37 GHz. The number of relatively high- 
luminosity HBLs increased with the higher assumed redshift. 
However, only at very high reshift values (z = 1 or z = 1.5) 
the luminosities of a few HBLs became roughly compara- 
ble to those of LBLs. Considering the population average, 
z av = 0.33, it is questionable whether all featureless sources 
would have z > 1, although for some of them this may be the 
case. Therefore, while the existence of high-luminosity HBLs 



in this sample is possible, we do not expect the correlations to 
be affected by them. 

In other frequency bands than radio, the effect of an in- 
creasing redshift was negligible. The shape of the correlations 
did not change notably, only the scatter increased somewhat. 
We note that in all wavelengths and for all assumed redshift 
values the significance of the statistical correlations remained 
the same. 



6. Other properties of the sample 

6.1. Log (S JS r ) -distribution 

The log (S X /S ^-distribution has traditionally been used to 
point out the bimodality in the BL Lac population. The divid- 
ing line between RBLs and XBLs has been log (S x /S r )=— 5.5 
jLaurent-Muehleisen et alJll999l) . when the fluxes are in the 
same units and X-ray and radio frequencies are 0.1-2.4 keV 
and 5 GHz respectively. Here we calculated the distribution to 
check the assumed continuity of the sample. X-ray fluxes are 
from the ROSAT band 0.1-2.4 keV and radio measurements 
are from Metsahovi at 37 GHz. If there was more than one flux 
measurement from one frequency, the average value was used. 

The distribution is indeed unbroken with an average of 
log (S x /Sr)=— 5.85 (fig. 01. The most populated interval is 
log (S x / S r )=[— 6, -6.5]. Pertaining to the unbalanced detection 
rates at 37 GHz, LBLs are overrepresented compared to HBLs 
in the plot. Only 14 % of the objects are HBLs. When con- 
sidered separately, LBLs, IBLs and HBLs move progressively 
from low to higher values of log{S X /S r ). There is substan- 
tial overlap between the distributions. The average values are 
-6.35, -5.50 and -4.85 for LBLs, IBLs and HBLs respec- 
tively. 
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Fig. 8. The colour plane: a ro plotted against a ox . Frequency in- Fig. 9. The colour plane: a ro plotted against a ox . Frequency in- 
tervals 5 GHz-5500 A-l keV, ROSAT band. Different symbols tervals 37 GHz-5500 A-l keV, ROSAT band. Different sym- 



denote the observational classification of BL Lacs. 



bols denote the observational classification of BL Lacs. 



6.2. Broad band spectral indices 
6.2.1 . The a ro - a ox diagram 

We calculated broad band spectral indices a ro and a ox to plot 
the sample in a ro vs. a ox diagram. Here the spectral index is 
defined as 



= - 



log(vi/v 2 ) 



(2) 



where S i and S 2 are the fluxes in frequencies v\ and V2 respec- 
tively and S v — v~". 

The a m vs. a ox diagram has also been used as a means to 
demonstrate the division of the population. Generally, XBLs 
have lower values of both indices occupying the lower left 
corn er of the plot, whereas RBLs lie in the upper rig ht cor- 
ner JStocke et alJl985ULaurent-Muehleisen et allll999l) . IBLs 
seem to have bridged the gap having intermediate values of 
both indices. 

For comparison, we plotted two diagrams; one with low 
radio frequency (5 GHz) (fig. [8} and one with high radio fre- 
quency (37 GHz) (fig. [5). The optical wavelength used was 
5500 A. In both diagrams we used X-ray data at 1 keV from 
literature and, in addition, data from ROSAT band 0. 1-2.4 keV. 
The error produced by the bandwidth difference is small com- 
pared with the benefits of a larger number of datapoints. 

When the lower radio frequency is used (fig.|S}, the indices 
for the whole population are a ox — 0.54-1.8 and a, o = 0.16- 
0.95. RBLs and XBLs occupy their expected locations but 
overall the distribution is even. 

As for the higher radio frequency diagram (fig.|9}, the inter- 
vals are slightly different; a ox = 0.63-1.8 and a ro — 0.29-1.13. 
The values of a ox are the same as in fig. |SJ but the number 
of datapoints is smaller. Here we note again the lack of HBLs 
compared to LBLs. This is why the lower left corner of fig. [9] 
is underpopulated compared with fig. [8] 

We also plotted the diagram with 5 GHz showing the 
LBL/IBL/HBL classification of our study (fig. [TO}. It clearly 
shows how the transition of the synchrotron peak from LBL 
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Fig. 10. The colour plane: a ro plotted against a ox . Frequency 
intervals 5 GHz-5500 A-l keV, ROSAT band. Different sym- 
bols denote the physical classification of BL Lacs obtained in 
this study. 

to HBL moves the object on the a ro - a ox -plane from the 
top to lower right and onw ards to lower left, as described by 
IPadovani & Giommll (fl995). However, there are several LBLs 
that appear on the wrong side of the a rx - 0.75 line that has 
usually been thought of as a dividing line between LBLs and 
HBLs. This suggests that the a rx = 0.75 divide is not very ef- 
fective and LBLs take on very scattered values of a ro . 

There are some objects in the a ro - a ox -plot (fig.llO> with 
no classification (marked with x). They are concentrated in 
quite a small area with a ox = 0.61—1.17 and a™ = 0.37-0.67. 
Their synchrotron peak frequencies have not been calculated 
on account of very poor fits, but judging by their spectral in- 
dices they are likely to be HBLs. 

6.2.2. a rx VS.\OQv peak 

The relation between a rx and log v pea i is shown in fig. ^2 
We used radio frequency 5 GHz and X-ray data comes pri- 
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Fig. 11. Radio-X-ray spectral index a rx plotted against syn- Fig. 12. Radio-optical spectral index a ro plotted against syn- 



chrotron peak frequency. 



chrotron peak frequency. 



marily from 1 keV. For those objects that have no 1 keV flux 
available we used data from the ROSAT band. According to 
the Spearman rank correlation test, the negative correlation is 
significant at 99 % level. However, the correlation seems to 
break apart at logy,,,^ < 14. LB Ls take on a rx -values be- 
tween 0.45-1. IFossati etal1 ( ll998h suggested that the correla- 
tion between the radio luminosity and X-ray luminosity in the 
rising Compton-component makes a rx tend to a fixed value 
when log v pea k < 14, but in our plot this seems not to be the 
case. 

When log Vp ea k > 14 the correlation is apparent. This speaks 
for the authenticity of the UHBL candidates, as they take their 
place at the end of the continuum with low spectral index val- 
ues. All objects with log v pea ^ > 19 also have a rx < 0.65. 

From the a rx vs.\ogv pea k -plot we can see that assigning 
an object an LBL/HBL classification based only on the value 
of a rx is risky. While object with a rx > 0.75 is very likely an 
LBL, the opposite does not hold. In fact, approximately 30 % 
of objects with a rx < 0.75 have log v pea k < 16. At a rx ~ 0.6 the 
possible values of v pea k span across as much as eight magni- 
tudes. 

When compar ed with the correspond ing figures of 
iFossati et alJ dl998b and lPadovani et alJ J2003I) . our plot seems 
to be a combination of the two. This is explained by the fact that 
both the samples used by Fossati et al. (Uy and Slew Survey) 
and part of the sample used by Padovani et al. (DXRBS) are in- 
cluded in the Metsahovi sample. Thereby we see a strong cor- 
relation with a substantial increase in scatter of the ff„- values 
towards the lower end. We stress that our data considerably ex- 
tends the range of values of log v pea k- 

6.2.3. a ro vs. log v peak 

The correlation between the spectral index between radio (5 
GHz) and optical (5500 A) frequencies and synchrotron peak 
frequency is shown in fig. El The overall trend is easily seen, 
changing from the steep negative correlation of LBLs to an al- 
most constant trend of HBLs. The change occurs roughly when 
the synchrotron peak of the SED moves to frequencies higher 



than those used to calculate a ro . After that point, the peak fre- 
quency no longer has an influence on a ro . The negative corre- 
lation is significant at 99 % level. Again, however, we see a few 
points breaking this scenario, both LBLs with very low values 
of a ro and HBLs with high values of a ro . 

We note that, as in the case of a rx , the scatter in the figure is 
significant. Therefore the value of a ro cannot be reliably used 
to classify BL Lacs. Only if a ro is very close to one, the object 
is likely to be an LBL. 

7. Discussion 

Our findings support the results of ICostamante et all J200ll) . 
Their study concludes that objects 1ES 0033+595, 1ES 
0120+340, 1ES 1218+304 and 1ES 1426+428 are ex- 
treme HBLs with synchrotron peak frequencies exceeding 
10 I8 Hz. Our results show that they all peak at roughly 
\ogv pea k ^19. M eanwhile, the extreme nature of 1ES 
2344+514 ( Giom mi et al l200ol) is not revealed in our study as 
it had a peak frequency of log v pea k = 16.4. Giommi et al. state 
that it is indeed very variable, with synchrotron peak frequen- 
cies ranging from log v pea i « 17 to log v pea i « 18. Mrk 501 is 
another example of a variable source; its synch rotron peak has 
been as high as log v pea t « 19 (Pian et al. 1998), but the peak in 
our SED occured at log v pea k ~ 16.84. 

In this paper we tested the correlation between the syn- 
chrotron peak frequency and luminosities at radio, optical and 
X-ray wavelengths and luminosity at the synchrotron peak of 
each sour ce. According to th e blazar sequence scenario pro- 
moted by IFossati et alJ dl998h . HBLs have lower luminosities 
at all wavelenghts. In the X-ray region, in the SED of LBLs the 
synchrotron and SSC-components meet resulting in a concave 
X-ray spectrum. In HBLs the synchrotron component peaks at 
or near the X-ray energies. Because the SEDs of LBLs and 
HBLs are so differently shaped in this frequency band, the 
X-ray correlation may be less significant. Still, on the whole, 
LBLs are expected to be more luminous. When the defined fre- 
quency bands are considered, our data seems to support this 
scenario in all other frequencies except in X-rays. There we 
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find a clear positive correlation which contradicts the blazar 
sequence. Luminosities rise towards the high peak frequencies. 
When all luminosity vs. \ogv pea k correlations are studied si- 
multaneously, a kind of sequence is revealed: at radio frequen- 
cies the negative correlation is steep, in optical wavelengths the 
correlation is only slightly negative and in the X-ray region it 
turns positive. However, the negative correlation of luminosity 
at 5 GHz and synchrotron peak frequency is expected to orig- 
inate solely from the shifting of the bulk of the synchrotron 
emission to the higher frequencies with growing log v pea k and 
as such does not authenticate the interdependece of source lu- 
minosity and peak frequency. A negative luminosity correlation 
at 5 GHz and a positive one at X-rays is what we would expect 
if the only changing parameter in the SED sequence was v pea k. 
When source luminosities at synchrotron peak frequencies are 
calculated, we find that they do not correlate with log v pea k. 

In addition to the positive correlation of X-ray luminosity 
and synchrotron peak and the lack of synchrotron luminosity 
correlation, the blazar sequence scenario is brought into ques- 
tion by the appearance of numerous low-luminosity LBLs at 
5 GHz. The faintest LBLs reach even lower luminosities than 
most of the HBLs. As the 5 GHz flux limit of even the least lim- 
ited radio-selected sample, RGB, is 20 mJy, it is obvious that 
these objects have been identified for the most part from X-ray 
surveys. We note that there are some 5 GHz low-luminosity 
LBLs that have a detection at 37 GHz. Especially for these 
sources we have to consider the possibility that the 5 GHz flux 
measurement is from a particularl y quiescent stat e and d oes not 
represent the object accurately. Pad ovani et alJ J2Q03I) found 
low-luminosity LBLs among the DXRBS BL Lacs. In our low- 
luminosity LBLs there are 4 objects from the DXRBS survey. 

The 5 GHz low-luminosity LBLs appear also in luminosity 
plots from other frequency bands, again, in the low-luminosity 
end. All the extreme HBLs have relatively high or intermediate 
luminosities. This effect, especially in X-rays, seems to sup- 
port the scenario that UHBLs are in fact a rare class. If they 
were common and X-ray luminous, the numerous X-ray sur- 
veys performed lately should have detected them in large num- 
bers. On the whole, the scatter in luminosities diminishes to- 
wards higher peak frequencies at all wavelengths. There is no 
solid evidence of truly high-luminosity HBLs with a measured 
redshift i n radio wa velengths. The high-luminosity HBLs re- 
ported bv lGiommi et al.l J2005t) are from the Sedentary Survey 
which is not included in our sample. 

We find that the low-luminosity LBLs behave abnor- 
mally also on the broad band spectral index plots, a ro vs. a ox , 
a rx \s.\ogv pea k and a ro vs. log v pea k- In the first the BL Lac 
classes take their expected places except for a few of the low- 
luminosity LBLs that invade the region usually populated by 
HBLs. Doing so, they cross the a, x = 0.75 border frequently 
used to classify BL Lacs to LBLs and HBLs. 

When the a rx vs. log v pea k -plot is examined, we find a 
strong negative correlation, but LBLs (log v pea k < 14) disrupt 
the trend and take on valu es a rx — 0.45-1 . Hints o f such a distri- 
bution was discovered bv lPadovani et alJ J2003h . but the wider 
range of log v pea k in our sample makes it evident. All LBLs with 
conspicuously low a rx are among the low-luminosity LBLs in 
the 5 GHz luminosity plot. We considered the possibility that 



the low indices in these objects are due to the wider bandwidth 
of the ROSAT fluxes compared to the 1 keV fluxes, but this 
seems to have little effect. For those objects that had a rx com- 
puted for both X-ray energies, we calculated the average val- 
ues of both indices. The average for ROSAT band a rx was 0.63 
while for 1 keV fluxes it was 0.69. This difference is small com- 
pared to the deviation of the low-a„ LBLs from their expected 
location. Also for 3 of the 10 LBLs with lowest a„-values 1 
keV flux has been used in calculating the spectral index. 

In a ro vs. log v pea k -plot LBLs behave in a similar way, 
at \ogv pea k ~ 13 they have a ro = 0.16-0.95. In this plot we 
also see 3 HBLs with conspicuously high a ro . They have rel- 
atively high radio luminosities and low optical luminosities. 
Each source has only one flux measurement at 5500 A, pos- 
sibly from an anomalously quiescent state. For these HBLs 
we also considered the possibility that their X-ray flux actu- 
ally originates from IC rather than synchrotron radiation. This 
would lower their peak frequencies and move them left on the 
a ro vs. log v pea k -plane. For two of them (RXS J1456. 0+5048, 
logv peflt =19.94 and RXS J1410.5+6100, log v peak =2Q.25) this 
seems possible albeit uncertain. 

As explained before, when there are few datapoints repre- 
senting the SED, the peak frequency is easily overestimated 
in the case of extreme HBLs. This affects also the luminosity 
correlations. However we do not expect the significance of the 
correlations at defined wavelengths to change, even if the data- 
points in the high-energy end moved to the left. The only case 
where the consequences could be substantial is the log L pea k vs. 
log v P eak -plot (fig.|6). If the peak frequency is exaggerated, so 
is the luminosity at peak frequency. Thus the uncertain data- 
points can move from upper right to lower left, possibly chang- 
ing the correlation. We tested this by artificially lowering the 
peak luminosities of those sources whose vF pea k was notably 
higher in the SED than their X-ray flux, vF x-ray We changed 
their peak frequencies accordingly, and the highest log v pea k in 
the sample changed from 21.46 to 19.44. Yet the overall shape 
of the plot did not change enough to produce significant corre- 
lation. Therefore we find it unlikely that the overestimation of 
the log v P eak of HBLs affects the main results of this study. 

8. Conclusions 

In this paper we collected a large amount of multifrequency 
data as well as new flux measurements at 37 GHz and plot- 
ted the spectral energy distributions of over 300 BL Lacs of 
the Metsahovi BL Lac sample. Using such an extensive sample 
allowed us to detect the whole range of synchrotron peak fre- 
quencies reaching up to the MeV-region. The main conclusions 
are as follows: 

1. For 22 objects we find that \ogv pea k > 19. This high- 
frequency tail of the v pea k distribution is unrepresented in most 
previous studies, but we find it essential for determining the 
properties of the population accurately. 

2. The positive correlation of X-ray luminosity and syn- 
chrotron peak frequency and the lack of correlation between the 
source luminosity at synchrotron peak and peak frequency con- 
tradict the blazar sequence scenario. The differences between 
LBLs and HBLs cannot be attributed to unequal luminosities. 
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3. When the broad band spectral indices a rx and a ro are 
plotted against log v pra j. we find substantial scatter in the fig- 
ures. This implies that the values of a rx and a ro cannot be reli- 
ably used in BL Lac classification. 

4. Based on the smooth, declining distribution of v pea k and 
the lack of any bimodality in all other tests and calculations, we 
conclude that the BL Lac population as a whole is continuous 
and undivided. 
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Table 3. The Metsiihovi BL Lac sample. The synchrotron peak frequency and subsequent classification designated in this study are included 
when available. 



Source 


KA.(JZUUU ) 


DeC.(JZUUU^ 


log V pea k 


/"'lace 

uass 


NRAO 5 


00:06:13.9 


-06:23:36 


12.75 


LBL 


RX J0007.9+4711 


00:07:59.9 


47:12:07 


16.14 


IBL 


MS 0011.7+0837 


00:14:19.7 


08:54:04 


16.74 


HBL 


RXS J00 18. 4+2947 


00:18:27.8 


29:47:32 


- 


- 


PKS 0017+200 


00:19:37.9 


20:21:46 


13.08 


LBL 


PKS 0019+058 


00:22:32.5 


06:08:05 


13.19 


LBL 


RXS J0325.2+1515 


00:35:14.9 


15:15:04 


13.73 


LBL 


1ES 0033+595 


00:35:52.6 


59:50:04 


18.93 


HBL 


1ES 0037+405 


00:40:13.8 


40:50:04 


16.8 


HBL 


RXS J0045.3+2127 


00:45:19.1 


21:27:43 


16.89 


HBL 


B3 0045+395 


00:47:55.2 


39:48:57 


- 


- 


EXO 0044.4+2001 


00:47:08 


20:17:44 


16.05 


IBL 


PKS 0047+023 


00:49:43.3 


02:37:04 


13.56 


LBL 


PKS 0048-097 


00:50:41.2 


-09:29:06 


13.39 


LBL 


NPM1G -09.0033 


00:56:20 


-09:36:32 


- 


- 


RXS J0058.2+1723 


00:58:16.8 


17:23:14 


- 


- 


QJ0109+181 


01:09:08.1 


18:16:03 


- 


- 


NPM1G +41.0022 


01:10:04.8 


41:49:50 


17.7 


HBL 


RXS J01 11.5+0536 


01:11:30.1 


05:36:28 


- 


- 


S2 0109+22 


01:12:05.8 


22:44:39 


13.59 


LBL 


RXS J01 15.7+2519 


01:15:46.5 


25:19:57 


13.43 


LBL 


1ES 0120+340 


01:23:08.7 


34:20:51 


18.32 


HBL 


MS 0122.1+0903 


01:24:44.5 


09:18:49 


15.53 


IBL 


B3 0133+388 


01:36:32.6 


39:06:00 


16.59 


HBL 


PKS 0139-09 


01:41:25.8 


-09:28:43 


13.4 


LBL 


1ES 0145+138 


01:48:29.7 


14:02:18 


15.44 


IBL 


NPM1G +01.0067 


01:52:39.6 


01:47:17 


- 


- 


8C 0149+710 


01:53:25.8 


71:15:07 


14.75 


IBL 


RXS J0155.9+1502 


01:56:00.3 


15:02:13 


- 


- 


87GB 01569+1032 


01:59:34.4 


10:47:07 


15.56 


IBL 


RXS J0200.4+2712 


02:00:29.5 


27:12:36 


- 


- 


MS 0158.5+0019 


02:01:06.1 


00:34:00 


17.87 


HBL 


RXS J0202.4+0849 


02:02:26.4 


08:49:12 


- 


- 


S5 0159+72 


02:03:33.3 


72:32:53 


13.99 


LBL 


MS 0205.7+3509 


02:08:38.2 


35:23:13 


15.22 


IBL 


87GB 02109+5130 


02:14:17.9 


51:44:52 


17.69 


HBL 


RXS J0216.5+2314 


02:16:32.1 


23:14:47 


- 


- 


Z 0214+083 


02:17:17 


08:37:03 


15.23 


IBL 


S4 0218+35 


02:21:05.4 


35:56:15 


- 


- 


3C 66A 


02:22:39.6 


43:02:08 


15.63 


IBL 


RXS J0227.2+0201 


02:27:16.6 


02:01:58 


- 


- 


1ES 0229+200 


02:32:48.6 


20:17:17 


19.45 


HBL 






lA'AI-^A 


1 S RA 


IRI 

LdL^ 


AO 0235 + 164 


02:38:38.8 


16:36:59 


13.57 


LBL 


S5 0238+71 


02:43:31 


71:20:18 


16.3 


IBL 


NPM1G +10.0097 


02:45:13.5 


10:47:23 






RXS J0250.6+1712 


02:50:38 


17:12:08 






MS 0257.9+3429 


03:01:03.8 


34:41:01 


13.28 


LBL 


4C 47.08 


03:03:35.2 


47:16:17 


14.18 


LBL 


RXS J0303.5+0554 


03:03:30.1 


05:54:17 






PKS 0306+102 


03:09:03.6 


10:29:16 


12.94 


LBL 


VZw331 


03:13:57.9 


41:15:24 


14.48 


LBL 


RXS J03 14.0+2445 


03:14:02.7 


24:44:31 


12.67 


LBL 



E. Nieppola et al.: SEDs of BL Lacertae Objects, Online Material p 3 



Table 3. continued. 



Source 


RA.(J2000) 


Dec.(J2000) 


log v peak 


Class 


RXS J03 14.3+0620 


03:14:23.9 


06:19:57 


19 


57 


HBL 


RXS J03 16. 1+0904 


03:16:12.9 


09:04:43 


15 


91 


IBL 


MS 03170+1834 


03:19:51.8 


18:45:35 


16 


99 


HBL 


RGB J032 1+2336 


03:22:00 


23:36:11 






- 


2E 0323+0214 


03:26:13.9 


02:25:14 


19 


87 


HBL 


RXS J033 1.3+0654 


03:31:19.4 


06:54:28 






- 


RXS J0349.9+0640 


03:49:59.7 


06:40:56 






- 


PKS 0406+121 


04:09:22.1 


12:17:39 


13 


22 


LBL 


2E 0414+0057 


04:16:52.4 


01:05:24 


20 


71 


HBL 


1WGA J0421.5+1433 


04:21:33.1 


14:33:54 


13 


93 


LBL 


MS 0419.3+1943 


04:22:18.5 


19:50:53 


16 


82 


HBL 


PKS 0420+022 


04:22:52.2 


02:19:27 






- 


PKS 0422+004 


04:24:46.8 


00:36:07 


15 


69 


IBL 


MCG 38364 


04:25:51.3 


-08:33:38 


17 


17 


HBL 


2EG J0432+2910 


04:33:37.7 


29:05:56 


14 


09 


LBL 


1ES 0446+449 


04:50:07.3 


45:03:12 






- 


PKS 0459+135 


05:02:33.2 


13:38:11 


13 


55 


LBL 


Q 0458+6530 


05:03:03.4 


65:34:10 


18 


12 


HBL 


RXS J0505.5+0416 


05:05:34.7 


04:15:54 


16 


94 


HBL 


1ES 0502+675 


05:07:56.1 


67:37:24 


19 


18 


HBL 


S5 0454+84 


05:08:42.5 


84:32:05 


13 


58 


LBL 


MG 0509+0541 


05:09:25.9 


05:41:35 


15 


34 


IBL 


4U 0506-03 


05:09:39 


-04:00:36 


17 


94 


HBL 


2E 0514+0626 


05:17:04 


06:29:39 






- 


1ES 0525+713 


05:31:41.7 


71:22:17 






- 


TEX 0554+534 


05:58:11.6 


53:28:19 


14 


44 


LBL 


MS 0607.9+7108 


06:13:42.8 


71:07:29 


14 


85 


IBL 


87GB 06216+4441 


06:25:18.3 


44:40:02 


13 


61 


LBL 


1ES 0647+250 


06:50:46.5 


25:03:00 


18 


28 


HBL 


B3 0651+428 


06:54:43.5 


42:47:59 


15 


12 


IBL 


NPM1G +42.0131 


06:56:10.6 


42:37:02 


17 


25 


HBL 


EXO 0706.1+5913 


07:10:30.1 


59:08:21 


21 


05 


HBL 


RXS J0712.3+5719 


07:12:18.7 


57:19:22 






- 


S5 0716+714 


07:21:53.3 


71:20:36 


14 


46 


LBL 


RXS J0723.2+5841 


07:23:13.2 


58:41:23 






- 


FIRST J0724.7+2621 


07:24:42.8 


26:21:30 


16 


39 


IBL 


PKS 0723-008 


07:25:50.7 


-00:54:56 






- 


FIRST J0730.4+3307 


07:30:26.1 


33:07:22 


16 


29 


IBL 


RXS J0737.3+3517 


07:37:21 


35:17:41 


17 


77 


HBL 


FIRST J0738.6+3 139 


07:38:37.8 


31:39:30 






- 


PKS 0735+17 


07:38:07.4 


17:42:19 


13 


95 


LBL 


FIRST J0741. 3+2253 


07:41:18.8 


22:53:39 






- 


MS 0737.9+7441 


07:44:05.1 


74:33:59 


13 


61 


LBL 


S4 0749+54 


07:53:01.3 


53:53:00 


13 


12 


LBL 


GB 0751+485 


07:54:45.7 


48:23:51 


14 


32 


LBL 


PKS 0754+100 


07:57:06.7 


09:56:35 


13 


63 


LBL 


RXS J0800. 1+6210 


08:00:06.5 


62:10:12 








RXS J0801. 0+6444 


08:01:00.7 


64:44:43 








RXS J0805.4+7534 


08:05:26.5 


75:34:25 


15 


96 


IBL 


SBS 0802+596 


08:06:25.9 


59:31:07 


16 


69 


HBL 


B2 0806+31 


08:09:13.4 


31:22:22 








RXS J0809.6+3455 


08:09:38.5 


34:55:37 


18 


29 


HBL 


1ES 0806+524 


08:09:49.2 


52:18:58 


16 


56 


HBL 


PKS 0808+019 


08:11:26.6 


01:46:52 


13 


17 


LBL 



E. Nieppola et al.: SEDs of BL Lacertae Objects, Online Material p 4 



Table 3. continued. 



Source 


RA.(J2000) 


Dec.(J2000) 


log V pe ak 


Class 


RXS J0812.1+5717 


08:12:08.8 


57:17:34 


- 


- 


EXO 0811.2+2949 


08:14:21.8 


29:40:32 


- 


- 


1WGA J08 16.0-0736 


08:16:04.3 


-07:35:57 


14.19 


LBL 


RXS J0816.3+5739 


08:16:23.8 


57:39:03 


17.19 


HBL 


RXS J08 16.6+6208 


08:16:40.9 


62:08:44 


15.25 


IBL 


OJ425 


08:18:16.1 


42:22:46 


13.33 


LBL 


FIRST J0818.4+2814 


08:18:27.3 


28:14:02 


16.01 


IBL 


FIRST J0819.4+4037 


08:19:25.8 


40:37:43 


16.61 


HBL 


FIRST J0820.3+3640 


08:20:20.2 


36:40:04 


14 


LBL 


4C 22.21 


08:23:24.8 


22:23:03 


13.09 


LBL 


PKS 0823+033 


08:25:50.3 


03:09:24 


13.08 


LBL 


RXS J0828.2+4153 


08:28:14.2 


41:53:50 


18.66 


HBL 


B3 0827+395 


08:30:19.4 


39:23:47 


- 


- 


PKS 0829+046 


08:31:48.9 


04:29:39 


13.53 


LBL 


1H 0827+089 


08:31:54.8 


08:47:58 


14.22 


LBL 


OJ448 


08:32:23.2 


49:13:21 


13.06 


LBL 


TEX 0836+182 


08:39:30.7 


18:02:47 


14.55 


IBL 


FIRST J0847.0+41 17 


08:47:02.5 


41:17:57 


18.11 


HBL 


RXS J0847.2+1133 


08:47:12.9 


11:33:52 


19.13 


HBL 


RXS J0848.4+8111 


08:48:27.8 


81:11:47 


- 


- 


US 1889 


08:54:09.9 


44:08:31 


17.42 


HBL 


OJ287 


08:54:48.8 


20:06:30 


13.89 


LBL 


NPM1G -09.0307 


09:08:02.2 


-09:59:37 


15.52 


IBL 


B2 0906+31 


09:09:53.3 


31:06:03 


17.4 


HBL 


Ton 1015 


09:10:37.1 


33:29:24 


15.39 


IBL 


FIRST J0910.8+3902 


09:10:52 


39:02:02 


18.33 


HBL 


B2 0912+29 


09:15:52.4 


29:33:24 


16 


IBL 


RXS J0916.8+5238 


09:16:52 


52:38:28 


17.22 


HBL 


MS 0922.9+7459 


09:28:02.6 


74:47:19 


- 


- 


RXS J0929.2+5013 


09:29:15.4 


50:13:35 


14.59 


IBL 


S5 0916+86 


09:29:42.7 


86:12:21 


14.16 


LBL 


1ES 0927+500 


09:30:37.5 


49:50:25 


21.13 


HBL 


B2 0927+35 


09:30:55.3 


35:03:38 


14.8 


IBL 


RXS J0930.9+3933 


09:30:56.8 


39:33:33 


- 


- 


SBS 0936+522 


09:39:37.9 


52:01:46 


- 


- 


B2 0937+26 


09:40:13.6 


26:03:26 


14.75 


IBL 


US 1015 


09:50:11.8 


45:53:20 


15.48 


IBL 


RGB J0952+656 


09:52:32.2 


65:38:01 


15.08 


IBL 


MS 0950.9+4929 


09:54:09.8 


49:15:00 


16.92 


HBL 


S4 0954+65 


09:58:47.2 


65:33:54 


13.76 


LBL 


4C 22.25 


10:00:21.36 


22:33:07.4 


- 


- 


MS 0958.9+2102 


10:01:42.4 


20:48:18 


15.54 


IBL 


EXO 1004.0+3509 


10:06:56.3 


34:54:45 


16.92 


HBL 


RXS J 1008 1+4705 


10:08:11.3 


47:05:20 


19.67 


HBL 


NRAO 350 


10:12:13.3 


06:30:57 


16.09 


IBL 


RXS J1012.7+4229 


10:12:44.3 


42:29:57 


20.97 


HBL 


GB 1011+496 


10:15:04.2 


49:26:01 


16.74 


HBL 


RXS J1016.2+4108 


10:16:16.7 


41:08:13 


16.62 


HBL 


RXS J1022.7-0112 


10:22:43.9 


-01:12:56 


17.97 


HBL 


1ES 1028+511 


10:31:18.5 


50:53:36 


18.56 


HBL 


FIRST J1032.6+3738 


10:32:40.7 


37:38:26 






RXS J1037. 7+5711 


10:37:44.2 


57:11:57 


14.95 


IBL 


TEX 1040+244 


10:43:09 


24:08:35 


13.1 


LBL 


1ES 1044+549 


10:47:45.8 


54:37:41 


13 


LBL 



E. Nieppola et al.: SEDs of BL Lacertae Objects, Online Material p 5 



Table 3. continued. 



Source 


K/Y.^JZUUU ) 




log Vp ea k 


i^i ass 


MS 1050.7+4946 


10:53:44.2 


49:29:54 


15.29 


IBL 


FIRST J1054.5+3855 


10:54:31.8 


38:55:22 


16.68 


HBL 


RXS J1055.5-0126 


10:55:34.1 


-01:26:05 


- 


- 


RXS J1056. 1+0252 


10:56:06.6 


02:52:13 


- 


- 


FIRST J1057. 3+2303 


10:57:23.1 


23:03:18 


18.7 


HBL 


RXS J1057.8+0059 


10:57:52.4 


00:59:13 


- 


- 


FIRST J1058.4+2817 


10:58:29.9 


28:17:46 


18.37 


HBL 


RXS J1058.6+5628 


10:58:37.7 


56:28:12 


15.64 


IBL 


MC 1057+100 


11:00:20.2 


09:49:35 


- 


- 


RXS Jl 100.3+4019 


11:00:21.1 


40:19:28 


18.76 


HBL 


RXS Jl 102.8-0148 


11:02:52 


-01:48:51 


- 


- 


MRK 421 


11:04:27.2 


38:12:32 


18.49 


HBL 


1ES 1106+244 


11:09:16.2 


24:11:20 


16.91 


HBL 


OP 1106.7+3654 


11:09:33.5 


36:38:26 


- 


- 


FIRST Jl 110.9+3539 


11:10:56.9 


35:39:06 


- 


- 


RXS Jl 110.6+7133 


11:10:37.5 


71:33:57 


16.96 


HBL 


RXS Jill 1.5+3452 


11:11:30.9 


34:52:00 


- 


- 


FIRST Jl 117.6+2548 


11:17:40.4 


25:48:46 


15.71 


IBL 


EXO 1118.0+4228 


11:20:48.1 


42:12:13 


17.41 


HBL 


US 2504 


11:29:50.1 


26:52:53 


12.97 


LBL 


MS 1133.7+1618 


11:36:17.6 


16:01:53 


15.9 


IBL 


MRK 180 


11:36:26.5 


70:09:28 


18.61 


HBL 


RXS Jl 136.5+6737 


11:36:30.1 


67:37:04 


17.55 


HBL 


FIRST Jl 136.8+2550 


11:36:50.1 


25:50:52 


15.13 


IBL 


2E 1 146+2456 


11:49:29.9 


24:38:55 


17.87 


HBL 


EXO 1449.9+2455 


11:49:30.3 


24:39:27 


19.83 


HBL 


B2 1147+245 


11:50:19.2 


24:17:54 


13.95 


LBL 


RXS Jl 151.4+5859 


11:51:24.6 


58:59:14 


16.4 


IBL 


FIRST Jl 152.1+2837 


11:52:10.7 


28:37:21 


14.8 


IBL 


FIRST Jl 153.7+3823 


11:53:42.9 


38:23:06 


- 


- 


MS 1154.1+4255 


11:56:46.6 


42:38:10 


14.91 


IBL 


B3 1159+450 


12:02:08.6 


44:44:21 


15.92 


IBL 


B3 1206+416 


12:09:22.8 


41:19:41 


14.59 


IBL 


1207+39W4 


12:10:26.6 


39:29:08 


- 


- 


MS 1209.0+3917 


12:11:34.2 


39:00:55 


- 


- 


1ES 1212+078 


12:15:10.9 


07:32:03 


15.91 


IBL 


Q 1214+1753 


12:16:56.9 


17:37:12 


- 


- 


B2 1215+30 


12:17:52 


30:07:01 


15.58 


IBL 


GB2 1217+348 


12:20:08.4 


34:31:22 


14.46 


LBL 


PG 1218+304 


12:21:21.9 


30:10:37 


19.14 


HBL 


ON 231 


12:21:31.7 


28:13:58 


14.84 


IBL 


UM493 


12:22:06.5 


-01:06:38 


- 


- 


RXS J1222.2+3541 


12:22:12.4 


35:41:00 


12.81 


LBL 


S5 1221+80 


12:23:40.4 


80:40:04 


14.21 


LBL 


MC 1991 C_i_9/l<9 
IVlo 1ZZ 1 . o-rZ^f jZ 


IZ.Z^f .Z^f . J 


Z^f . jO.Z^f 


1 ^ QQ 


T RI 


1WGAJ1225. 3+1818 


12:25:18.2 


18:18:20 


14.71 


IBL 


FIRST J1226.0+2604 


12:26:04.1 


26:04:28 


16.29 


IBL 


RXS J1230.2+2517 


12:30:14 


25:18:06 


14.9 


IBL 


2E 1258+1437 


12:31:23.9 


14:21:25 


14.91 


IBL 


MS 1229.2+6430 


12:31:31.3 


64:14:18 


16.25 


IBL 


B2 1229+29 


12:31:43.6 


28:47:49 






FIRST J1236.3+3900 


12:36:23.1 


39:00:01 


16.61 


HBL 


MS 1235.4+6315 


12:37:38.6 


62:58:44 


15.98 


IBL 


RXS J1237. 0+3020 


12:37:05.5 


30:20:04 







E. Nieppola et al.: SEDs of BL Lacertae Objects, Online Material p 6 



Table 3. continued. 



Source 


RA.(J2000) 


Dec.(J2000) 


log V pe ak 


Class 


1ES 1239+069 


12:41:48.3 


06 


36 


01 


17.38 


HBL 


RXS J1241.6+3440 


12:41:41.2 


34 


40 


32 


- 


- 


Ton 116 


12:43:12.7 


36 


27 


44 


- 


- 


PG 1246+586 


12:48:18.8 


58 


20 


29 


15.05 


IBL 


1ES 1249+ 174E 


12:51:45.5 


17 


11 


17 


- 


- 


FIRST J1252.3+2640 


12:52:19.5 


26 


40 


53 


- 


- 


S4 1250+53 


12:53:11.9 


53 


01 


11 


14.82 


IBL 


1ES 1255+244 


12:57:31.9 


24 


12 


40 


16.89 


HBL 


MS 1256.3+0151 


12:58:54.6 


01 


34 


43 


- 


- 


MS 1258.4+6401 


13:00:17.6 


63 


44 


39 


16.35 


IBL 


FIRST J1301.7+4056 


13:01:45.7 


40 


56 


24 


16.55 


HBL 


RXS J 1302.9+5056 


13:02:55.5 


50 


56 


17 


- 


- 


MC2 1307+12 


13:09:33.9 


11 


54 


24 


13.07 


LBL 


1WGAJ1309.6+0828 


13:09:38.9 


08 


28 


28 


14.64 


IBL 


AUCVn 


13:10:28.66 


32:30:43.8 


13.75 


LBL 


HS 1309+2605 


13:12:19.2 


25 


49 


58 


- 


- 


TEX 1312+240 


13:14:43.8 


23 


48 


26 


15.84 


IBL 


RXS J1319.5+1405 


13:19:31.7 


14 


05 


34 


20.85 


HBL 


UM566 


13:19:55.1 


01 


52 


58 


- 


- 


1ES 1320+084N 


13:22:54.9 


08 


10 


10 


13.12 


LBL 


RXS J1324.0+5739 


13:24:00.8 


57 


39 


16 


15.48 


IBL 


RXS J1326.2+2933 


13:26:15 


29 


33 


29 


- 


- 


RXS J1326.2+1230 


13:26:17.6 


12 


30 


00 


16.32 


IBL 


RX J 1340. 1+2743 


13:40:10.9 


27 


43 


48 


- 


- 


RXS J 1340.4+44 10 


13:40:29.8 


44 


10 


04 


16.51 


HBL 


RXS J1341.0+3959 


13:41:05 


39 


59 


45 


20.06 


HBL 


1338.8+2705 


13:41:05.8 


26 


50 


26 


- 


- 


1340.8+2721 


13:43:05.1 


27 


06 


24 


- 


- 


RXS J1353.4+5601 


13:53:28 


56 


00 


55 


19.23 


HBL 


FIRST J 1354.4+3706 


13:54:26.7 


37 


06 


54 


16.92 


HBL 


RXS J1359. 8+5911 


13:59:53.7 


59 


11 


01 


13.66 


LBL 


MC 1400+162 


14:02:44.5 


15 


59 


57 


16.47 


IBL 


RXS J1404.8+6554 


14:04:49.6 


65 


54 


30 


- 


- 


MS 1402.3+0416 


14:04:51 


04 


02 


02 


15.83 


IBL 


MS 1407.9+5954 


14:09:23.5 


59 


39 


41 


16.63 


HBL 


PKS 1407+022 


14:10:04.6 


02 


03 


07 


13.69 


LBL 


RXS J1410.5+6100 


14:10:31.7 


61 


00 


10 


20.25 


HBL 


FIRST J1414. 1+3430 


14:14:09.3 


34 


30 


57 


- 


- 


RGB J1415+485 


14:15:36.8 


48 


30 


30 


14.56 


IBL 


PKS 1413+135 


14:15:58.8 


13 


20 


24 


12.83 


LBL 


CRSS 1416.3+1137 


14:16:20.7 


11 


37 


37 


- 


- 


2E 1415+2557 


14:17:56.6 


25 


43 


25 


19.24 


HBL 


OQ530 


14:19:46.6 


54 


23 


14 


14.16 


LBL 


RXS J1422.6+5801 


14:22:39 


58 


01 


55 


- 


- 


FIRST J1426. 1+3404 


14:26:07.7 


34 


04 


26 


14.1 


LBL 


PKS 1424+240 


14:27:00.5 


23 


48 


00 


15.7 


IBL 


RGB J1427+541 


14:27:30.3 


54 


09 


24 


14.89 


IBL 


H 1426+428 


14:28:32.7 


42 


40 


21 


18.55 


HBL 


TEX 1428+370 


14:30:40.6 


36 


49 


03 


14.26 


LBL 


CSO 474 


14:36:45.7 


35 


57 


01 






RXS J1436.9+5639 


14:36:57.8 


56 


39 


25 


17.5 


HBL 


PG 1437+398 


14:39:17.5 


39 


32 


43 


16.7 


HBL 


1ES 1440+122 


14:42:48.3 


12 


00 


40 


16.45 


IBL 


MS 1443.5+6349 


14:44:34.9 


63 


36 


06 


17.05 


HBL 
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Table 3. continued. 



Source 


RA.(J2000) 


Dec 


(J2000) 


log V pe ak 


Class 


RXS J1445.0-0326 


14:45:06.2 


-03 


:26 


:12 


- 


- 


RXS J1448.0+3608 


14:48:00.6 


36 


08 


31 


16.73 


HBL 


RXS J1449.5+2746 


14:49:32.8 


27 


46 


21 


18.82 


HBL 


RXS J1451 4+6354 


14:51:27.5 


63 


54 


19 


- 


- 


SBS 1452+516 


14:54:27.1 


51 


24 


33 


- 


- 


RXS J1456.0+5048 


14:56:03.7 


50 


48 


25 


19.94 


HBL 


RXS J1458.4+4832 


14:58:28 


48 


32 


40 


21.46 


HBL 


B3 1456+375 


14:58:44.8 


37 


20 


22 


13.47 


LBL 


MS 1458.8+2249 


15:01:01.9 


22 


38 


06 


15.26 


IBL 


FIRST J1502. 1+2528 


15:02:08.3 


25 


28 


45 


14.17 


LBL 


FIRST J1502.5+3350 


15:02:34 


33 


50 


55 


13.66 


LBL 


RXS J1508. 7+2709 


15:08:42.9 


27 


09 


10 


17.3 


HBL 


SBS 1508+561 


15:09:48 


55 


56 


17 


15.22 


IBL 


FIRST J15 15.9+2426 


15:15:56.2 


24 


26 


20 


15.89 


IBL 


RXS J1516.7+2918 


15:16:41.6 


29 


18 


10 


18.62 


HBL 


PKS 1514+197 


15:16:56.8 


19 


32 


12 


13.6 


LBL 


1H 1515+660 


15:17:47.6 


65 


25 


24 


18.11 


HBL 


FIRST J1530.7+2310 


15:30:44 


23 


10 


13 


- 


- 


FAQS J1530.7+5329 


15:30:44.5 


53 


29 


28 


16.88 


HBL 


RXS J1532. 0+3016 


15:32:02.2 


30 


16 


29 


17.05 


HBL 


RXS J1533.4+3416 


15:33:24.3 


34 


16 


41 


18.32 


HBL 


RGB J1534+372 


15:34:47.2 


37 


15 


55 


14.26 


LBL 


1ES 1533+535 


15:35:00.8 


53 


20 


37 


19.68 


HBL 


FIRST J1535.4+3922 


15:35:29.1 


39 


22 


46 


16.38 


IBL 


MS 1534.2+0148 


15:36:46.8 


01 


37 


59 


18.83 


HBL 


1ES 1544+820 


15:40:15.8 


81 


55 


07 


17.79 


HBL 


4C 14.6 


15:40:46.5 


14:47:45.9 


14.85 


IBL 


RXS J1542.9+6129 


15:42:56.9 


61 


29 


56 


14.72 


IBL 


RXS J1544.3+0458 


15:44:18.7 


04 


58 


24 


16.77 


HBL 


MS 1552.1+2020 


15:54:24.1 


20 


11 


25 


17.12 


HBL 


PG 1553 + 11 


15:55:43.1 


11 


11 


24 


16.49 


IBL 


MYC 1557+566 


15:58:48.5 


56 


25 


14 


- 


- 


RXS J1602.2+3050 


16:02:18 


30 


51 


09 


16.42 


IBL 


PKS 1604+159 


16:07:06.4 


15 


51 


34 


14.73 


IBL 


RXS J1610.0+6710 


16:10:04.1 


67 


10 


26 


17.45 


HBL 


B3 1615+412 


16:17:06.6 


41 


06 


45 


14.41 


LBL 


87GB 16166+2206 


16:18:47.9 


21 


59 


26 


15.66 


IBL 


4C 37.46 


16:21:11.3 


37:46:04.9 


- 


- 


RXJ 1624.7+3726 


16:24:43.4 


37 


26 


42 


- 


- 


RXS J1624.9+7554 


16:24:56.5 


75 


54 


55 


13.3 


LBL 


RXS J1626.4+3513 


16:26:25.6 


35 


13 


38 


15.28 


IBL 


RXS J1631. 3+4217 


16:31:24.7 


42 


17 


03 


18.99 


HBL 


RXS J1638.0+7326 


16:38:02.6 


73 


26 


10 


- 


- 


RXS J1644.2+4546 


16:44:20 


45 


46 


44 


17.48 


HBL 


1643.2+4021 


16:44:53.2 


40 


16 


28 






FIRST J1645.9+2947 


16:45:57.7 


29 


47 


30 






FIRST J165 1.1 +4212 


16:51:09.2 


42 


12 


53 


17.82 


HBL 


RXS J 165 1.6+72 18 


16:51:41.5 


72 


18 


19 






RGB J 1652+403 


16:52:50 


40 


23 


10 


14.97 


IBL 


MRK501 


16:53:52.2 


39 


45 


36 


16.84 


HBL 


B3 1659+399 


17:01:24.6 


39 


54 


36 


13.41 


LBL 


FIRSTJ1702.1+2643 


17:02:09.6 


26 


43 


15 


14.27 


LBL 


FIRST J1702.6+3 115 


17:02:38.6 


31 


15 


43 


13.49 


LBL 


RXS J1704.8+7138 


17:04:46.9 


71 


38 


18 


15.58 


IBL 
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Table 3. continued. 



Source 


RA.(J2000) 


Dec.(J2000) 


log V pe ak 


Class 


MS 1704.9+6046 


17:05:34.9 


60:42:17 


17.17 


HBL 


FIRST J1712.8+2931 


17:12:48.8 


29:31:17 


18.79 


HBL 


RXS J1718.6+7358 


17:18:40.1 


73:58:15 


- 


- 


PKS 1717+177 


17:19:13.1 


17:45:07 


13.08 


LBL 


RXS J1719.6+7443 


17:19:41 


74:43:00 


14.19 


LBL 


B2 1722+40 


17:24:05.5 


40:04:38 


13 


LBL 


H 1722+119 


17:25:04.4 


11:52:16 


15.8 


IBL 


IZwl87 


17:28:18.6 


50:13:11 


17.42 


HBL 


RXS J1728.6+7041 


17:28:38.3 


70:41:08 


13.6 


LBL 


RXS J1732. 0+6926 


17:32:05.4 


69:26:16 


- 


- 


OT465 


17:39:57.1 


47:37:59 


13.91 


LBL 


RGBJ 1742+597 


17:42:32 


59:45:06 


14.18 


LBL 


NPM1G +19.0510 


17:43:57.9 


19:35:09 


17.91 


HBL 


B3 1743+398B 


17:45:37.7 


39:51:31 


17.71 


HBL 


1ES 1745+504 


17:46:32.3 


50:28:09 


- 


- 


B3 1746+470 


17:47:26.6 


46:58:51 


13.96 


LBL 


S4 1749+70 


17:48:33.1 


70:05:50 


14.55 


IBL 


B3 1747+433 


17:49:00.4 


43:21:52 


13.56 


LBL 


RXS J1750.0+4700 


17:50:05 


47:00:44 


18.1 


HBL 


PKS 1749+096 


17 


51:32.7 


09:39:01 


12.78 


LBL 


RXS J1756. 2+5522 


17 


56:15.9 


55:22:18 


19.9 


HBL 


MS 1757.7+7034 


17 


57:13.3 


70:33:37 


13.7 


LBL 


S5 1803+784 


18 


00:45.4 


78:28:04 


14.05 


LBL 


RXJ1801.7+6638 


18 


01:46.7 


66:38:40 


- 


- 


3C 371 


18 


06:50.7 


69:49:28 


14.65 


IBL 


RGB J1808+468 


18 


08:01.2 


46:49:41 


14.63 


IBL 


RGB J181 1+442 


18 


11:53.5 


44:16:29 


15.07 


IBL 


B2 1811+31 


18 


13:35.3 


31:44:17 


15.53 


IBL 


4C 56.27 


18:24:07.07 


56:51:01.5 


12.95 


LBL 


RXS J1829.4+5402 


18 


29:24.3 


54:03:00 


15.07 


IBL 


Q 1832+6845 


18 


32:36.6 


68:48:09 


- 


- 


RXS J1838.7+4802 


18 


38:49.2 


48:02:34 


13.52 


LBL 


RGB J1841+591 


18 


41:20.3 


59:06:08 


14.91 


IBL 


RXS J1848.7+4245 


18 


48:47.1 


42:45:39 


- 


- 


1ES 1853+671 


18 


53:52.1 


67:13:55 


16.64 


HBL 


87GB 19021+5536 


19 


03:11.6 


55:40:39 


14.5 


LBL 


S4 1926+61 


19 


27:30.4 


61:17:32 


13.44 


LBL 


RXS J 193 1.1 +0937 


19 


31:09.5 


09:37:13 


- 


- 


1ES 1959+650 


19 


59:59.9 


65:08:55 


18.03 


HBL 


S5 2007+77 


20 


05:31.1 


77:52:43 


13.15 


LBL 


S5 2010+72 


20 


09:52.3 


72:29:19 


13.64 


LBL 


PKS 2012-017 


20 


15:15.1 


-01:37:33 


14.68 


IBL 


S5 2023+76 


20 


22:35.6 


76:11:26 


14.1 


LBL 


PKS 2032+107 


20 


35:22.3 


10:56:06 


14 


LBL 


1ES 2037+521 


20 


39:23.5 


52:19:50 


16.13 


IBL 


PKS 2047+039 


20 


50:06.2 


04:07:49 


13.83 


LBL 


S5 2051+74 


20 


51:33.8 


74:41:40 


18.75 


HBL 


PKS 2131-021 


21 


34:10.2 


-01:53:17 


12.76 


LBL 


MS 2143.4+0704 


21 


45:52.3 


07:19:27 


13.92 


LBL 


PKS 2149+17 


21 


52:24.8 


17:34:39 


13.85 


LBL 


BL LAC 


22 


02:43.3 


42:16:39 


14.28 


LBL 


RXS J2209.3+1031 


22 


09:18.5 


10:31:43 


13.36 


LBL 


RXS J2219.7+2120 


22 


19:45.3 


21:20:48 


17.61 


HBL 


RXS J2225.1 + 1136 


22 


25:11.2 


11:36:01 
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Table 3. continued. 



Source 






log Vp ea k 


i^i ass 


PKS 2223-114 


22:25:43.7 


-11:13:41 


14.65 


IBL 


3C446 


22:25:45.1 


-04:56:34 


13.37 


LBL 


S5 2229+69 


22:30:35.6 


69:46:29 


13.09 


LBL 


RXS J2233.0+1335 


22:33:00.9 


13:35:59 


16.61 


HBL 


RGB J2243+203 


22:43:54.7 


20:21:04 


14.15 


LBL 


B3 2247+381 


22:50:05.8 


38:24:37 


15.61 


IBL 


PKS 2254+074 


22:57:17.3 


07:43:12 


14.18 


LBL 


RXS J2304.6+3705 


23:04:36.6 


37:05:08 


21.01 


HBL 


B3 2311+396A 


23:13:50.2 


40:03:03 


- 


- 


QJ2319+161 


23:19:43.4 


16:11:51 


15.48 


IBL 


TEX 2320+343 


23:22:44 


34:36:14 


16.73 


HBL 


1FS 2321 +419 


23-23-54 1 


42:11:19 


13.26 


LBL 


B3 2322+396 


23:25:17.9 


39:57:37 


15.29 


IBL 


1ES 2326+174 


23:29:03.3 


17:43:30 


18.07 


HBL 


QJ2338+212 


23:38:56.4 


21:24:41 


17.62 


HBL 


MS 2336.5+0517 


23:39:07 


05:34:36 


14.91 


IBL 


1ES 2344+514 


23:47:04.8 


51:42:18 


16.4 


IBL 


MS 2347.4+1924 


23:50:01.7 


19:41:52 


15.8 


IBL 


RXS J2350.3-059 


23:50:17.9 


-05:59:28 






TEX 2348+360 


23:50:36.7 


36:22:11 


16.1 


IBL 


PKS 2354-02 


23:57:25.1 


-01:52:15 







E. Nieppola et al.: SEDs of BL Lacertae Objects, Online Material p 10 



-8 


- NRAO 5 




10 






12 




— 


14 






16 













-10 












-12 








-14 


> 








-16 



- RXJ0007.9+4711 




-8 


- MS 0011.7+0837 


i i i i 


^ - 


-2 - 

erg cm s 


-10 
-12 




• 

— 


,4 


> 

2 


-14 
-16 







10 12 14 16 18 20 22 
log V [Hz] 



10 12 14 16 18 20 22 
log V [Hz] 



10 12 14 16 18 20 22 
log V [Hz] 




10 12 14 16 18 20 22 
log V [Hz] 



10 12 14 16 18 20 22 
log V [Hz] 



10 12 14 16 18 20 22 
log V [Hz] 













-10 










bt) 


-12 








-14 


> 


M 






-16 



- EXO 0044.4+2001 




-8 






-10 




'a 






rg c 


-12 








: y 


> 


-14 






-16 




8 10 12 14 16 18 20 22 
log V [Hz] 



8 10 12 14 16 18 20 22 
log V [Hz] 



8 10 12 14 16 18 20 22 
logv [Hz] 



- NPM1G +41.0022 


1 








1 




10 12 14 16 18 20 22 
log V [Hz] 



10 12 14 16 18 20 22 
logv [Hz] 



8 10 12 14 16 18 20 22 
logv [Hz] 




- MS 0122.1+0903 








»' 





10 12 14 16 18 20 22 
log V [Hz] 



10 12 14 16 18 20 22 
logv [Hz] 




10 12 14 16 18 20 22 
logv [Hz] 



Fig. 13. Spectral energy distributions of the Metsahovi BL Lac sample. Only datapoints used in the fit are shown in the figure. 
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. 13. continued. 
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